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Abstract 



Next-to-next-to-leading order (NNLO) QCD corrections to Higgs boson 
hadroproduction have recently been calculated in the heavy top-quark limit 
m t — > oo. The m t —>■ oo limit introduces double-logarithmic corrections in 
Inx, with x = m 2 H /s. We identify these corrections order by order in a s . As 
an application, we derive an analytic expression for the dominant i C 1 part 
of the NNLO coefficient. 
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Recently, Harlander and Kilgore [|IJ have completed the calculation of 
the next-to-next-to-leading order (NNLO) QCD correction to Higgs boson 
production in hadron-hadron collisions. The calculation is done in the heavy 
top quark limit, in which m t — > oo and all other quark masses vanish. In 
this limit, the effective Lagrangian of is used for coupling the Higgs boson 
to gluons. The hard scattering coefficient expanded through NNLO has the 
structure 



C{a S) x,m 2 H /^ 2 ) = c (°)(x) + — c (1) (x) + c (1) (x) L 

+ [ C ( 2 \x) + c^\x) L + 5 (2) (x) L 2 ] + ••■ (I) 



where L = \n(m 2 H / fi 2 ) , with m# the Higgs boson mass and /i the factorization 
scale; x = mj^/s, with s the center-of-mass energy. The term in is 
the leading-order term J|; the term in S l \c^ is the next-to- leading-order 
term [|]. 

Ref. |1| presents the result for as an expansion in powers of (I — 
x), x — 1 being the production threshold. Explicit numerical results are 
given up to a very high order in this expansion, (1 — x) 16 . Once the hard 
scattering coefficient is convoluted with the parton densities, this expansion 
leads to a corresponding expansion for the cross section. This is demonstrated 
numerically to converge well [|IJ. As discussed by Catani, de Florian and 
Grazzini || and by Harlander and Kilgore the reason for the convergence 
of the (1 — x) expansion is mainly kinematical, and depends on the steeply 
rising behavior of parton luminosities for decreasing momentum fractions. 
This behavior also explains why results based only on soft || [7|] and virtual [[8] 
contributions give a reasonably good approximation to the full answer. 

In this note we focus on the influence of the mt —>■ oo approximation 
on the i C 1 part of the coefficient. It was observed in || that a local 
Higgs-gluon coupling gives rise to double- logarithmic corrections in lnx. In 
this note we identify these corrections explicitly order by order in a s . In 
particular we obtain an analytic result for the dominant x C 1 part of the 
NNLO coefficient. 

Although the i <C 1 behavior is unlikely to affect the convergence of 
the inclusive cross section at NNLO, it may however affect distributions as- 
sociated to less inclusive observables. This analysis may thus be useful for 
further evaluations of the hard scattering. 
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Let us consider the a; < 1 corrections in the framework of the effective 
Lagrangian [fj] for the Higgs boson coupling to gluons: 

£ e ff = (G F V2) 1/2 K(a s ) G°^ v G aiJ,u H , (2) 

where K(a s ) = a s /(12n) + 0(a 2 ) is the coefficient function containing the 
dependence on the top quark mass, and is known to order a A s |7]. The 
x C 1 corrections come, in higher order graphs, from integrating two-gluon 
irreducible amplitudes over the transverse momenta ki and k 2 of the glu- 
ons that couple to the Higgs boson. In the full theory, contributions from 
|ki|, |k 2 | m t are suppressed by the top quark loop. In the effective theory, 
large values of | ki | and | k 2 1 are allowed all the way up to the kinematic limit 
y^s. This results [[| in the perturbative coefficients having stronger singu- 
larities in the complex iV plane as iV — > (with iV the moment variable of 
Eq. (H)). More precisely, double poles a s /N 2 appear order by order in a s . 

To determine these contributions, recall the structure of the hard scat- 
tering coefficient in the MS subtraction scheme for iV — > 0: 

C^{a s ,m 2 H / fi 2 ) = / dxx N ~ l C{a s ,x,m 2 H / fi 2 ) 
Jo 

= R 2 N ( lN ) (m 2 H /n 2 ) 2 ^ h N ( lN , lN ) . (3) 
Here Rm is a normalization factor associated with the choice of the MS 



scheme [[LI]; 7^ is a known, universal anomalous dimension, with the per- 



turbation expansion 

the function is constructed from the matrix element a for the off-shell 
amplitude g(ki) + g{k2) — * H by taking the following integral transform 



^(71,72) = 7i 72 



d 2 k, ( k 2 \ 71 fd 2 k 2 ( k| V 2 



7rkf \m 2 H J J ixk.2 \m 2 H 







r kl 


k 2 \ 


X 


X, 






Jo X 




m H J 



(5) 

The matrix element a, although off-shell, is defined gauge-invariantly by cou- 
pling the gluons with eikonal polarizations, as in p. It is readily calculated 
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using the effective Lagrangian (0): 



a x 



ki k 2 \ a 2 s m 2 H G F \/2 m% (k x • k 2 ) 2 



' m H 'm H J 288tt x 2 k?k|. ( m % + kf + k 2 ) 2 



x yx m H / 

Now consider the function Hn- The x integral in Eq. (|5[) can be done 
using the 5-function of Eq. (|j): 

1 dx „ a 2 s m 2 H G F y/2 m% (k x • k 2 ) 2 
o IT X 288^ kfkf { m \ + kf + k 2 ) 2 



2 \iV-2 
H) 



\m 2 +(k 1 + k 2 ) 2 } N - 2 



(7) 



Notice that when iV — > the matrix element is unsuppressed for large 
|ki|, |k 2 |. This is a consequence of approximating the Higgs coupling by 
a local vertex, and is in contrast with the general case of a finite top-quark 
mass |L2|. The behavior at large transverse momenta in Eq. (^) causes the 
function h to have N = singularities not only through its dependence on 
the anomalous dimension @ but also through its explicit N dependence. 
Introducing the notation 

^ = \4/m 2 H , & = kj/m 2 H (8) 

and performing the angular integral in Eq. (^) we get 

- 2 - 2 G F V2 



h N {jhj2) = ^ 7i 72 / d£i£i / d& 6 

2887T Jo Jo 



72-1 
2 



2(i+a+6r v 2 2 ' '(i + 6 + e2) 2 y 

G6(iV-l)(iV-2) /iV iV+l 466 M (9) 

2(i + 6 + 6) 2+w 2 X V2' 2 ' d '(i + 6 + 6)VJ ' 1 J 

where 2 -F\ is the hypergeometric function. The integration region of large 6 
and 6 i n Eq. (§) gives rise to a pole at 71 + 7 2 = N. This pole is the origin 
of the double-logarithmic terms. 
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To approximate h near N = taking the contribution of this pole into 
account, we expand the hypergeometric functions in Eq. (Q) to the lowest 
order around N = 0, 



N-l N-2 



(JV-l)(JV-2) 2 F 1 



466 



66 



;i+6+6) 2 , 

'at iV + 1 „ 



1 + 6 + 6 



+ ©(iV) , (10) 



466 



;i + 6 + 6 



2 + 0(N) , (11) 



but retain the A" dependence in the factor (1 + 6 + 6) N ■ The integrations 
in 6 an d 6 can now be carried out simply in terms of Euler F functions: 



^(7i,72) 



288^^ 
1 



x 



r(i + 7l )r(i + 72 )r(iv-7 1 -7 2 ) 
1 



T(N) 



+ 37172 



T(N + 2) 



iV<l . 



(12) 



We can now expand h for small anomalous dimensions, setting 7l = 72 = 
7 at, with 7 tv given in Eq. (|), and determine the leading A^ 1 part of the 
perturbative coefficients through the second order: 



hN(lN,lN] 



a. s m H Gp 
2887^ 



(13) 



1 

TV 



1 + 27^+47^— 



+ 



■) ( 7JV < 1 , N < 1) 



Via Eq. @, Eq. (0) gives the double-pole terms a s /N 2 order by order in 
perturbation theory. These terms are generated by the pole (N — 2 7 tv) _1 in 
h, see Eq. ([TJ). 

We are now in a position to find the NNLO correction to the gluon-gluon 
coefficient for 3; < 1. Higher order contributions to K(a s ) [1(J in Eq. (0) 
are subleading for 1 < 1. The factor in Eq. (§) is a single- logarithmic 



factor |l (J that starts to contribute one order higher than NNLO: 

R N = 1 + 0(a s /Nf . (14) 
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The (m 2 H / n 2 ) 2lN factor in Eq. ([$]) is to be expanded for small 7 at, as has been 
done with h. Then, by taking an inverse transform from moment space to 
x-spa.ce, we get 



C gg (a s ,x,m 2 H / \i 2 



+ 



ot 2 m 2 H GF 
2887TV/2 

2 



r 2 



5(1 -x) + — C A (-2 \nx + 2 L) (15) 

7T 



In 3 x + 2 In 2 x L - 2 In x L 2 + 



The first two terms in the square bracket reproduce respectively the leading- 
order result H and the i < 1 piece of the next-to- leading-order result ||; 
the third term determines the i<l pieces of the coefficients c g 2 g \ c g 2 g \ c g 2 g 
in Eq. CUD- 

The quark components of the coefficient also contain double-logarithmic 
corrections. These can be determined in an analogous manner. We find 



C qg (a s ,x,m 2 H / '/j 2 ) 



ot 2 m 2 H GF 
2887rv/2 

3 



7T 



\nx + L)+(^j 2 C F C A 



-In x -\ — in x L hi x L 

2 2 2 



(16) 



C q q(a s ,x,m H /n ) 



x 



a 2 m 2 H GF 

®s\ 2 r2 



(17) 



• - In 3 x + In 2 x L 

3 



In a; L + 



In summary, perturbative corrections to Higgs production are double- 
logarithmic for x 1 in the m t — > oo effective theory, while single-logarithmic 
in the full theory. We have given simple formulas that incorporate all the 
double logarithms to any order in a s in the m t — > oo case. The formulas also 
contain subleading iV — > terms, mostly of kinematical origin, which may 
be relevant for matching the x < 1 contributions with expansions around 
the soft limit x — > 1 |], ||, p. For observables sensitive to the x <C 1 region, 
it will be relevant to include finite top-mass corrections. 
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